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We investigate a formulation of continuum 4d gravity in terms of a constrained BF theory, in the 
spirit of the Plebanski formulation, but involving only linear constraints, of the type used recently 
in the spin foam approach to quantum gravity. We identify both the continuum version of the linear 
simplicity constraints used in the quantum discrete context and a linear version of the quadratic vol- 
ume constraints that are necessary to complete the reduction from the topological theory to gravity. 
We illustrate and discuss also the discrete counterpart of the same continuum linear constraints. 
Moreover, we show under which additional conditions the discrete volume constraints follow from 
the simplicity constraints, thus playing the role of secondary constraints. 
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I. INTRODUCTION 

The equations of general relativity can be derived from several different action principles [lj], leading 
to equivalent classical theories (in the case of pure gravity, at least). Among them we can mention, in 
addition to the Einstein-Hilbert action Q , the Palatini first order formulation and its modification proposed 
by Hoist 0. This last one is of special interest because it is the classical, covariant starting point for the 
canonical quantisation leading to Loop Quantum Gravity 0, Even if not the only possible useful one 
[6] , a particularly popular action in covariant approaches to quantising gravity 0] , like the s pin foam Q 
and group field theory approach Q, is the formulation as a constrained BF (or Plebanski 10, 11 1) theory. 



Here one starts from topological BF theory [12j in 4 spacetime dimensions, and adds suitable constraints 
on the two-form B variables of the theory such that, on solutions of these constraints, the action reduces 
to the Palatini or Hoist action for general relativity. We will summarise the idea behind this formulation in 
the following. In the original Plebanski formulation the constraints on the B variables are quadratic, and 
so are the discrete constraints that are then implemented in the spin foam models based on a simplicial 
discretisation. On the other hand, the most recent developments in the spin foam and group field theory 
approach to quantum gravity are based on a linear set of discrete constraints, which can be shown to be 
slightly stronger, in the restrictions they impose on the original BF configurations, than some of the original 
discrete quadratic constraints. Once more, we will detail this construction in the following. 



In this note we investigate whether a formulation in terms of linear constraints is also possible in the 
classical continuum theory, and what it implies. We will see that the replacement of diagonal and cross- 
simplicity constraints with linear constraints at the continuum level is relatively straightforward, after one 
has introduced new variables ha forming a basis of three- forms at each point. One then needs additional 
constraints corresponding to the volume constraints. We will see that one can linearise these constraints 



* sg452@cam.ac.uk 

t daniele.oriti@aei.mpg.de 



2 



too. We then give a discrete version of these linear volume constraints, which bears a striking resemblance 
to the so-called "edge simplicity" constraints of 13J . We note that only certain linear combinations of the 
volume constraints one would naively write down are necessary to constrain the bivectors Y> AB (A) sufficiently. 
Similarly to the quadratic case, we will also see that when linear diagonal and off-diagonal constraints hold 
everywhere in a 4-simplex, and one also imposes closure constraints on both bivectors T, AB (A) and normals 
71.4(A) referring to the tetrahedra of the 4-simplex, the sufficient set of linear combinations of the linear 
volume constraints follows. This additional four-dimensional closure constraint on the normal vectors has, 
to the best of our knowledge, not been considered or implemented as an additional condition in the spin 
foam literature yet, although it does appear in some first order formulation of Regge calculus (l4| . and it 
plays also a role in the discrete analysis of 15 1. 



II. GRAVITY AS CONSTRAINED BF THEORY: CONTINUUM AND DISCRETE RESULTS 

Let us briefly review what is known at the classical continuum and discrete level, concerning the Plebanski 
formulation of classical gravity. We limit our considerations to the covariant, Lagrangian context, and to 
a very small subset of the available results, those which have been already of direct relevance for quantum 
gravity model building, especially in the s pin foam context. For recent results in the canonical Hamiltonian 
setting, both continuum and discrete, see 



Consider first the Einstein- Hilbert-Palatini- Hoist 0,0] Lagrangian (without cosmological constant), 1 

S = ITn I [\ ^abcd E A A E B A R cd [uj] + — E A A E B A R AB [u>]\ , (2.1) 

where spacetime is assumed to be of the form TZ x R so that a (3+1) splitting can be performed, u> AB is a 
G-connection one-form (the gauge group G is SO(3, 1) or SO (4), or an appropriate cover), R AB its curvature, 
and E A is an M 4 -valued one-form representing an orthonormal frame. The term involving 7, known as the 
Hoist term, is not relevant classically; its variation vanishes if torsion is zero (but see 19]). It is, however, of 
fundamental importance in loop quantum gravity (LQG), and more generally for any canonical formulation 
of gravity, as it modifies the symplectic structure of the theory. 
If one introduces a g-valued two-form, 

tAB _ 1 ( 1 AB TpC a rpD 1 1 rpA A 77-B 



8^gU 15 7 /' (2 ' 2) 

then the action (|2.1I) becomes 

S = f B AB A Rab M + X a C a [B] , (2.3) 



i.e. it takes the form of a topological BF theory with additional constraints C a which enforce that B AB is 
indeed of the form (12.2[) . and that are enforced by means of Lagrange multipliers A Q . 

As said, BF theory without constraints is topological. Its equations of motion imply that uj ab is flat and 
the covariant exterior derivative of B AB vanishes. Having no local degrees of freedom, the quantisation of 
such a theory is therefore rather simple and quite well understood. Inspired by this classical formulation, 
the main issue when trying to construct a quantum theory related to quantum gravity, in 4 dimensions 2 , is 



1 We follow the usual conventions: Capital Latin indices are internal indices, small Latin indices are coordinate indices. For 
an abstract collection of indices, we use small Greek letters. 

2 Notice that, in 3 spacetime dimensions, gravity in 1st order form coincides, classically, with 3d BF theory, and that in higher 
dimensions a similar formulation of gravity as constrained BF theory can be given 12011 . 
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then the correct implementation of appropriate constraints that lead to (|2.2[) for some set of one- forms E A , 
either at the level ofquantum states or in a path integral formulation. Indeed, the bulk of the work in the 
spin foam approach [8, 21 - 24 1 (as well as in the group field theory formalism 0, [2|| 26|), in recent years, has 



been devoted to this task. These constraints are also the subject of this note. 
To simplify the following calculations, we introduce another two- form field T, AB , 



^ab _ 



1 ( r>AB T AB rCd\ io a \ 

t^{ b ~r cdB )' (2 - 4) 



where s is the spacetime signature, s = — 1 for G = SO(3, 1) and s = +1 for G = 50(4) (and we assume 
7 2 7^ s) 3 . This is a linear redefinition which simplifies the constraint (|2.2I) . 



S AB = -^—E A AE B , (2.5) 

87T7G 

but leads to more terms in the action. The translation of all calculations from one set of variables to the 
other is usually straightforward. 

The traditional way to enforce the restriction (|2.5p . the one matching the original classical Plebanski 



formulation of gravity, was to add quadratic simplicity constraints 4 to the action 27], 1 2 81] , 

eABCD^ab^cF = Vtabcd, (2.6) 

where V can be expressed in terms of T, AB by contracting (I2.6[) with e abcd , to give: V = 
^ e abcd eABCD j]AB-£CD _ This is itself a reformulation of the original Plebanski constraint, which would 
read: 

abcd^ABsr^CD _ rr ABCD /r, 7 \ 

e ^ab ^cd v C , (Z.I ) 

and is equivalent to the first under assumption that V ^ everywhere. The version (|2.6p has the advantage 
of permitting a much simpler discretisation and thus a more straightforward implementation within the spin 
foam formalism. Under the same assumption V ^ 0, there are the following four classes of solutions to (|2.6j) : 

either Y, AB = ±e A A e B or Z AB = ±^e AB CD e c A e D (2.8) 

for some set of one- forms e A (the factor 8irjG can obviously be introduced by rescaling). One would like 
to select only the first class of solutions E AB = +e A A e B , which, when substituted in the BF action, gives 
the Hoist action (I2.1j) . Classically, this is not a severe problem. As shown in 27], non-degenerate initial 
data of a solution of the form T, AB = +e A A e B generically remain within the same branch of solutions. The 
situation in the quantum theory, where one necessarily has contributions from all branches, is less clear. 

More troublesome, if V = 0, the field Y, AB does not permit a straightforward geometric interpretation 
at all. Since in the region of the phase space where V = 0, the theory is less constrained, and hence has 
more degrees of freedom, these non-geometric configurations should be expected to be dominating in a path 
integral [271 ] . unless measure factors are such that this is avoided. 

Spin foam models are usually defined in a piecewise flat context, and spin foam amplitudes are defined 
for given simplicial complexes [8j. Therefore one is interested in identifying a discrete version of the above 
constraints that could be imposed at the level of each complex. The version (|2.7[) of the simplicity constraints 
admits only a rather involved discrete counterpart |28| and, upon quantisation, leads to the Reisenberger 



model [28|, |29| , which has so far received only limited attention. 



3 For uniformity of the discussion, we shall in the following talk about "time" and use the label even when the gauge group 
is SO (A) and the signature Riemannian. 

4 "Simplicity" because a two-form that can be written as a wedge product of one-forms is called simple. 
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The discrete analogue of the constraints (|2.6p led instead (271. |28| to the construction of the Barrett-Crane 
model [3o| . in the case in which the Immirzi parameter is excluded from the original action (7 — > 00). The 
construction is initially limited to a single 4-simplex, the convex hull of 5 points in R 4 (IR 1,3 , in the Lorentzian 
case) with the topology of a 4-ball, whose boundary is triangulated by the 5 tetrahedra identified by the 5 
independent subsets of 4 such points, while subsets of 3 points identify the four triangles belonging to each 
of these 5 tetrahedra, each of the triangles being shared by a pair of tetrahedra. One then associates a Lie 
algebra element (bivector) G so (4) ~ A 2 M 4 (similarly in the Lorentzian case) to each triangle A in a 
given triangulation by integrating the two-form H AB over A. The task is then to constrain appropriately 
these Lie algebra variables (or their quantum counterpart) following the continuum treatment. 

It is useful to split the set of continuum equations (12.61) into two sets. Out of the 21 equations (12.61) . one 
first identifies and imposes those 18 which have zero on the right-hand side (the "diagonal" and "off-diagonal" 
simplicity constraints) , 

eABCD^tb^ab 1 = e ABCD^ab^F = V a, 6, C. (2.9) 

This corresponds to the case if one or two of the indices of the two fields £ coincide. At the discrete level, 
this translates into two triangles on which the same fields are discretised which either coincide or at least 
share a single edge, and thus belong to the same tetrahedron. Thus all bivectors are required to satisfy 

^abcd^a^a = (diagonal simplicity constraint) 

and 

^abcd ^zy^ = for all A, A' sharing an edge (cross-simplicity constraint). 

These two sets of equations can be imposed at the level of each tetrahedron in the 4-simplex. 

The remaining three equations (the "volume" constraints) are equivalent to the requirement that: 

^ABCD^Ql^i^ = e ABCD^02^'li > = e ABCD^(yF ^23* ^ ^(^0> (2-10) 

and can be imposed at the discrete level as the requirement that, for each 4-simplex: 

£abcd ^zy 5 — y for all A, A' not sharing an edge (volume constraints) (2-11) 

where V is defined by the above equation, and is interpreted, on the solutions of the constraints, as the 
volume of the 4-simplex. 

An additional condition on the bivectors is usually considered, namely the "closure" constraint, which 
states that the sum of four bivectors corresponding to the faces of one tetrahedron is zero: 



= . (2.12) 
This constraint can be understood in two ways. One can either view it as the condition that the triangles 



AC* 



described by the variables close to form a tetrahedron 3l|, or as a consequence of the equations of 

motion. In a topologically trivial region such as the interior of a tetrahedron, a flat connection can be set 
to zero by a gauge transformation. Then using Stokes' theorem, the integral over the equation dY> AB = 
can be written as J A Y> AB = 0, which is the closure constraint. The canonical counterpart of this condition 
is then the Gauss constraint, which generates local gauge (rotation) transformations and is to be imposed 
on the quantum states of the theory. 

The same picture appears in three spacetime dimensions, where there are no simplicity constraints and 
one directly deals with a su(2) connection one- form oj a and an su(2)-valued one- form e A . Here the equation 
de A — is integrated over a (spacetime) triangle to give a closure constraint. The vectors (using su(2) ~ R 3 ) 
associated to the edges of the triangle add up to zero, and thus have a consistent geometric interpretation 
as edge vectors in R 3 . In this sense, an n-form with vanishing exterior derivative and appropriate internal 



5 



indices can be given a geometric interpretation as describing n-simplices closing up to form an (n+l)-simplex. 
We shall encounter another instance of this statement later on. 

The closure constraint, being linear in the £'s and local in each tetrahedron, is obviously easier to impose 
at the discrete level, and in the quantum theory, than the volume constraints. Thus it is a useful fact that it 
can indeed be imposed instead of them. More precisely, it can be shown 23j that the volume constraints in 
each 4-simplex are implied if one has enforced the diagonal and cross-diagonal simplicity constraints, plus the 
closure conditions everywhere, i.e. in all the tetrahedra of the 4-simplex (in general, i.e. for non-degenerate 
4-simplices, involving tetrahedra belonging to different "time slices"). From a canonical perspective, this 
observation is usually phrased as an interpretation of the volume constraints as "secondary constraints" 
required to guarantee conservation of the other constraints (including the Gauss (closure) constraint) under 
time evolution. 

After a period of investigations, several potentially worrying issues have been put forward regarding the 
Barrett-Crane model 32h34| (for a more recent analysis of the geometry of the Barrett-Crane model, see 



26[), and have given impetus to the development of alternative spin foam models 2l|, |22|, |24j . These models 
are known to have nice semiclassical properties 35|, and, importantly, generalise the spin foam setting to 
include the Immirzi parameter at the quantum level (for an early attempt, see [Ifjj]), and thanks to this allow 
for a more direct contact with the canonical loop quantum gravity. Their study is still somewhat preliminary, 
but the above properties make them promising candidates for a quantum theory related to gravity. One of 
the central features of the new models is the replacement of the quadratic simplicity constraints (|2.9[) by 
linear constraints of the form 



n y i(A)E" 4B (A) = VAC 



(2-13) 



where is the normal associated to the tetrahedron A and A is any of the faces of A. 

It can be shown that these are lightly stronger than the discrete diagonal and off-diagonal quadratic 
simplicity constraints, and remove some of the discrete ambiguity in the solution for £ : out of the classes 
of solutions (|2.8p . one can restrict to (the discrete version of) T, AB = ±e A A e B only. For a geometric 
analysis of these conditions in the discrete setting, see lj| 21, 24, 37]], and for a proof that the same discrete 
conditions can also lead to the Barrett-Crane model, see [2 



III. LINEAR CONSTRAINTS FOR BF-PLEBANSKI THEORY 



The purpose of this note, as anticipated, is to investigate whether a formulation in terms of linear con- 
straints is also possible in the classical continuum theory, and what it implies. 

Let us work backwards, at first. Assume that the two-form field T, AB is of the form T, AB — e A A e B , and 
that the "frame field" e A is non-degenerate, i.e. that the matrix (e A ) is invertible. It follows that 

e%E AB = S c a e B -S c b e B . (3.1) 

In order to make a connection to the discrete setting it is more convenient to work with exterior powers of 
the cotangent bundle only (n- forms can be integrated over n-dimensional submanifolds) . Hence we multiply 
(|3.1[) by e c def and insert the relation e ct je/e^ = (det e a A ) £ADEFe d e E ej , which is true for invertible matrices, 
obtaining 

e A DEF e£ef ef Y, AB = (det e A ) (e adef e B - e bdef e B ) . (3.2) 

One can define the three-form nAdef = n-Aldef] = ^ADEFe d e E e F , so that (|3.2j) take the form 

n Adef Y, A b B = (det e A ) (e arfe /ef - e 6de/ ef ) . (3.3) 

nAdef can be interpreted as a 3d volume form for the submanifold parametrised by (x d , x e ,x^) embedded in 
4d spacetime, whose internal index gives the normal to this submanifold. If e A are a basis of one-forms at 



G 



each spacetime point, then ua are a basis of three-forms at each spacetime point, and so one can choose to 
work either with one or the other set of variables. Clearly e A can be reconstructed from ua'- 



^e cdef n Ad ef = s(dete A )e c A = s^/det ^-e bd -fn Bdef j e% (3.4) 

This means that the set of variables n A (x) define a co-tetrad frame at any point of the spacetime manifold 
(for the discrete analogue of the above, see fisl]). 



A. Linearised Diagonal and Off-Diagonal Constraints 

So far we have just rewritten the equation we want to obtain for T, AB . Let us now consider the implications 
of imposing ()3.3j) as constraints, where we restrict to those with zero right-hand side, i.e. those for which 
{a, b} C {d, e, /}. These are half of the equations (|3.3[) . This will identify the continuum analogue of the 
linear simplicity constraints. 

Claim 1. For a basis ua of three- forms, the general solution to 

n Ad ef^ab B =0 V{a,b}c{d,e,f} (3.5) 

is 

Z* b B = G ab e?e?, (3.6) 

where e A is defined in terms of riAdef as in (|3. 41) . and assumed to be non-degenerate, and G a b = Gba 
and G aa = 0. Obviously, as the variables E and the tetrad field e A , the "coefficients" G a b are spacetime 
dependent. 

Proof. First note that we can rewrite (|3.5[) as 

D EF^AB n /q t\ 

£ADEFe d e e e f l, ab = (3.7) 
with e A defined by (I3.4[) . Then e A by assumption defines a basis in the cotangent space, so that 

^b B = G 9 a h b e A e B (3.8) 
for some coefficients G 9 ^ with G 9 ^ = G^^. Substituting this into (|3.7[) . we get 

^ e ADEF e A e%e*e F f G%e% = e gdef det(e A )e% G%, (3.9) 
and since det(e^) ^ and ef form a basis of (the internal) R 4 , this implies that 

Me/G^ = V{a,&}c{d,e,/}. (3.10) 
It follows that G 9 ^ = unless {g, h} — {a, b} and so G 9 ^ = 5^ 9 d^G a b. 

□ 

By a linear redefinition e A = X a e A one might try to set some of the G a b to a given value (usually ±1, but 
one might prefer ± 8?r G ) , but it is clear that one needs two conditions on the G a b for this to be possible. 

In the discrete context, one sets n^i(A) = (1,0,0,0) for each A by a gauge transformation. 5 One could 
use some of the gauge freedom here to restrict the form of ha' This amounts to finding a convenient 



5 This presumably involves an implicit assumption, namely that there is a non-degenerate normal to each tetrahedron, as well. 
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parametrisation for the coset space GL(A)/SO(3,l). Let us make the (usual) assumption that the normal 
to hypersurfaces {t = constant} is indeed timelike. Then one can use the boost part of 5*0(3, 1) to set 
^Ai23 = (0,0,0,0). The remaining 50(3) subgroup can then be used to make the (3 x 3) matrix n/ 0£ ; e , 
where I £ {1,2,3}, upper diagonal, so that one has the form 



/* 


* 


* 







* 


* 


* 








* 


* 













(3.11) 



Clearly, when this form of riAdef is assumed, integrating the three-form ua over a region where t is constant 
gives a vector in IR that only has a time component. Its magnitude specifies the three-dimensional volume 
of such a region. 



B. Linearised Volume Constraints 

As in the quadratic case, further constraints, in addition to the linear simplicity constraints (|3.5p . are 
needed to complete the identification T, AB = ±E A A E B . 
First of all, one can show the following. 

Claim 2. Under the assumption that all G a b are non-zero, the necessary and sufficient conditions for the 
existence of a linear redefinition E A = X a e A , such that either T, AB = cE A A E B or Y> AB = -cE A A E B , 
where c is a given positive number, are 

G12G03 = G01G23 = Gi2,Go2{^ 0). (3.12) 

Proof. Set c = 1; the extension to arbitrary c amounts to a further rescaling by \fc. Then the required 
redefinition is possible if and only if there exist Aq, . . . , A3, such that either G a b — A a At> for all a 7^ b or 
Gab = — A a Ab for all a 7^ b. Clearly (13. 12)) are necessary. They are also sufficient: Take 



Ai 



G12G13 



G 2 



. Gi2Gi3\ G12 . / Gi2Gi3\ G13 

A 2 = sgn(— — , A 3 = sgn — — , (3.13) 

(^23 J M \ ^23 / -^1 



which solves the equations for G12, G13 and G2 3 with sgn ^ G g,^ 13 specifying the overall sign. The remaining 
three equations for Gqi, G02 and G03 are then solved by the two relations (|3.12p and 

/ G12G13 \ G i , . 

A = sgn^— j— . (3.14) 

□ 

The assumption G a b 7^ is necessary: One solution to p.!2j) is Gi 2 = G23 = G13 = with the other G a b 
non-zero, which cannot be expressed as G a b — ±A a Ah. 

Further constraints, in addition to the linear simplicity constraints (|3.5j) . are needed to complete the 
identification Y> AB — ±E A A E B . One possibility is to use the quadratic volume constraints p,12[) . Take 
the three volume constraints (|2.10|) . 

vABriCD _ ^AB^CD , ^AB^CD /o K i 

CABCD^oi ^23 — tABCD^m ^13 = ^ABCD^oi ^23 1 [6. 10) 

and substitute the solution T, AB = Gabelfe^ of (|3.5|) . This gives precisely (|3.12|) . The non-degeneracy 
assumption needed for (j3. 12[) is then the usual one, namely V ^ in (|2.6p . 

This shows that imposing the linear version of the diagonal and off-diagonal simplicity constraints (|3.5|) 
together with the quadratic volume constraints (|2.10[) and a non-degeneracy assumption on Yi AB implies 
that 

S AB = ±cE A A E B (3.16) 
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for some set of one- forms E A , where c > can be chosen at will. Thus, linearising the diagonal and off- 
diagonal simplicity constraints means that two of the four types of solutions for Y> AB are removed, but on 
the other hand one needs to introduce a basis of three- forms ua at each spacetime point, which is put in as 
an additional variable. One also still has to assume 7^0. 

There is also generically no evolution of initial data with V ^= into a degenerate Y> AB with V — and a 



non-geometric interpretation (this is part of the discussion of 271]). The geometry of the spacetime manifold 
is specified by E A and not by e A which is only used to determine normals in the constraints. 

Alternatively, one might prefer to use a linear version of the volume constraints as well. Consider the 
original equation 



UAdef^ab = (dete a ) (e adef e b - e bdef e a ) , (3.17) 

which was equivalent to T, AB = e A A e B for an invertible frame field. So far we only considered one half of 
these equations, namely those with {a, 6} C {d, e, /}. The other half have the form 

nAbef^ab = (detef)e a i, e /ef , no sum over b, (3.18) 

with tabef 7^ 0. One way to read these equations is as the requirement on the left-hand side to be totally 
antisymmetric in (a, e, /): 

nAbef^ab B = n Abfa^tb B = n Abae^f B ■ (3.19) 

We could again try to turn the argument around and impose (|3.19[) as constraints on a g- valued two-form 
Y* AB together with the linear simplicity constraints (|3.5p . Substituting the solution ^ AB = G a b^ae^} of the 
linear simplicity constraints into p.!9[) gives (after diving by a non-zero factor ^ det(e„ )) 

£af)e/G Q bef = £ebJ a G e be B — 6fbaeG fb^b ■ (3.20) 

For e a bef 7^ this would imply G a b = G e b = G/&. By Claim 2, imposing (|3.19l) for one fixed 6, say b = 0, is 
generically not sufficient: If we know that Goi = G02 = G03 7^ 0, we still have the condition 

G12 = G13 = G23, (3.21) 

so that one would need more conditions of the form (|3.19[) . These will then imply that all G a b are equal, 
Gab = c 1 for some c' that could be positive, negative, or zero. One can absorb |c'| by an overall redefinition, 
so that one has 

YsAB =±cE A AE B , (3.22) 

for any chosen c, as before. Note that here it is possible, if d — at a point, that all E A are zero this point 
and so T, AB = as well. While this is a very degenerate geometry, it is still a geometry. 

While the conditions p,19[) . imposed for all values of b, are therefore sufficient to complete the identification 
of the two-form field T, AB as ±cE A A E B , note that (|3.19|) is a massively redundant set of constraints: In 
order to obtain at most five relations on the coefficents G a b (two relations (|3.12[) if all G a b are nonzero), we 
are imposing eight vector equations! We have not exploited the fact that (|3.20[) is a multiple of one of the 
vectors ef , which are by assumption linearly independent. We could add several of the equations (|3.19j) for 
different 6, instead of considering all equations for different b separately. Let us try to impose 

n Ab ef^ab B = 0, e €{0,1,2, 3} fixed. (3.23) 

b {a,f}<Z{b,e} 

Again substituting the solution T, AB — G a fceL A e^' of the linear simplicity constraints into f|3 . 23[) . we obtain 

~det(e A )J2 ^b e f G ab e B =0, e fixed, (3.24) 

6 {a,m{b,e} 
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which implies, by linear independence of the e B , that indeed G a t, = Gfb for all e ^ {a, b, /}. It is then 
sufficient to impose the constraint (|3.23[) for three different choices of e, say e = 0, 1, 2, so that we only need 
three vector equations instead of eight. 

By absorbing the constant G a t, — d (all G a b are equal) we rescale all e A by the same factor to obtain 
the variables E A that will have the physical intepretation of frame fields encoding the metric geometry of 
spacetime. While in the case of quadratic volume constraints the one-forms e A , or alternatively the three- 
forms n A , only specified the normals to submanifolds {x a — constant}, for linear volume constraints they 
can be directly interpreted, up to a position-dependent normalisation, as specifying an orthonormal basis in 
the cotangent space. 

Note that this implies that one can assume a convenient normalisation for the one- forms e A . Instead of 
just assuming non-degeneracy det(e^) ^ 0, one could fix dct(e^) = 1. This is no restriction of the physical 
content of the theory as the e A , for both linear and quadratic volume constraints, only have a geometric 
interpretation after rescaling. One could then interpret e A as a map into SX(4,R). For linear volume 
constraints, the relation between the normalised one-forms e A and the variables E A that are interpreted as 
frame fields is a single function on spacetime which may be viewed as a "gauge" in the sense of Weyl [38| . 

In contrast to the case of the quadratic volume constraint, no non-degeneracy assumption on the two- 
form Y, AB is needed to enforce simplicity. One might get T, AB = in some region as a solution to the 
constraints, in which case the action for this region will be zero. This is analogous to a metric with vanishing 
determinant in general relativity and, in contrast to the requirement V ^ outlined above, not an additional 
issue. Notice, however, that one still has to assume that the tetrad field e A and, equivalently, the co-tetrad 
field n A are non-degenerate, in order for the simplicity and volume constraints to imply p.!6[) . Failing this, 
one gets solutions of the constraints that admit no proper geometric interpretation. 



In the end, writing the action for BF theory in terms of £ 



AB 



S j B AB A Rab — / ^ AB A Rab + \t AB CD YF u A R AB , (3.25) 



we substitute (|3.16[) into this action, which gives (setting c = g^ 1 G ) 

S=^ f <t(x) (\e ABCD E A A E B A R CD + -E A A E B A R Al 



(3.26) 



One is left with a field cr(x) that can take the values ±1, but in the classical theory one may again argue 
that if a — 1 everywhere on an initial hypersurface, there will be no evolution into a = — 1. What we obtain 
is first order general relativity where one uses (dete) = ±| dete| instead of | dete| as a volume element in 
the action. If a is continuous as classical fields usually are assumed to be, this differs from the action with 
| det e| by an overall sign at most. 

To summarise, we have identified both a linear version of the quadratic simplicity constraints and a linear 
version of the (quadratic) volume constraints in the continuum, which can be used to reduce topological 
BF theory to 4d gravity in the continuum. We have found also that both linear versions are slightly 
stronger (i.e. more restrictive) than the corresponding quadratic constraints, so that the resulting constrained 
theory is likely to be closer to gravity at the quantum level than the one in which quadratic constraints are 
implemented. We now discuss the discrete counterpart of the constraints found above. 

C. Discrete Linear Constraints and their Relations 



The discrete analogue of (|3.5[) is just the linear constraint used in [2l|, l24( , as desired: 

n A (A)Z AB {A) = VAcA. (3.27) 
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One could write down also a discrete version of (|3.19p , obtained in the natural way, demanding that within 
the same 4-simplex 

n A (A)Z AB {A') = n A {A')H AB (A") (3.28) 
whenever A' <£_ A and A" <£_ A' and the edge shared by A' and A is the same as that shared by A' and A". 



In the following we adopt the notation of [23j, where the tetrahedra in a given 4-simplex are labelled by 
A, B, C, D, E, so that triangles are represented by AB, AC, etc., and edges by combinations ABC, ABD, 
etc. The orientation of the triangles and tetrahedra in (I3.28|) is then fixed by the signs of the permutations 
of the letters, 

n A (A A )E AB (A BC ) = -n A (A B )S AB (A AC ) = n A (A c )Z AB (A AB ), etc. (3.29) 

In analogy to the continuum case, it will be sufficient to impose, instead of the full set of conditions (|3.29p . 
certain linear combinations of p.29[) to complete the geometric interpretation of the bivectors T, AB (A). 
The discrete analogue of the three continuum equations (|3.23|) . where the index e was kept fixed, is to pick 
one of the tetrahedra and add those six equations out of (I3.29j) which involve triangles belonging to this 
tetrahedron. Starting with A, we impose the constraint 

J2 nA(Ai)S AS (A Aj ) = 0, (3.30) 

and the equivalent conditions for the tetrahedra B to E, thereby needing to satisfy only five instead of 20 
volume constraints. 

The above discrete formulation of the linearised volume constraints resembles strongly the edge simplicity 
constraints studied, in a canonical setting, in (li| . and it imposes indeed the same restriction on the discrete 
data. However, it does not match exactly any of the various expressions given for these edge simplicity 
constraints in [13]. The correspondence between the two, therefore, deserves to be studied in more detail, 
given also that edge simplicity constraints have been shown to be crucial for the kinematical phase space of 
BF theory (and of loop gravity) to reduce to that of discrete gravity, in accordance with what we find here 
in a covariant setting. 



In spin foam models such as [2l|,[24|, as we mentioned earlier, only the diagonal and off-diagonal simplicity 
constraints, but no quadratic volume constraints (|2 . 10[) are imposed. This is because in the discrete setting, 
one can use the closure constraint (I2.12p , imposed in all the tetrahedra in a 4-simplex, to relate the (quadratic) 
simplicity constraints to the volume constraints, so that if the former are imposed everywhere the latter 
follows. Since the quadratic simplicity constraints follow from the linear ones, as can be easily checked, this 
argument is still valid if one uses linear simplicity constraints. 

One might hope that the sufficient set of linear volume constraints (|3.30[) would also follow from the linear 
simplicity constraints and the closure constraints. This is almost the case, but not quite. In fact, one more 
constraint should be added to simplicity and closure imposed in the five tetrahedra in the 4-simplex. This 
is a "4d closure" constraint of the form 

n A {A A ) + n A {A B ) + n A (A c ) + n A (A D ) + n A {A E ) = 0, (3.31) 

where A, are the (appropriately oriented) tetrahedra of a given 4-simplex. 

Just as for the usual closure constraint (|2.12[) . there are two ways to understand why such a constraint 
must be imposed. Recall that if one demands the triangles described by discrete variables S^ 8 close to form 
a tetrahedron, they have to satisfy (|2.12[) . Alternatively, one can start with the continuum field equation 
Vr^E^ = 0, where V' w ' is the covariant derivative for the connection oj ab , set the (flat) connection to 
zero by a gauge transformation, and integrate this over an infinitesimal 3-ball (whose triangulation is a 
tetrahedron). 

The new constraint (|3.31|) seems to be the analogous statement that tetrahedra close up to form a 4- 
simplex. By Hodge duality A 1 ^ 4 ~ A 3 R 4 and any internal covector n A (A) can be mapped to a three-form; 
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unlike for two- forms, any three- form can be written as e 1 Ae 2 Ae 3 for some e a . Demanding that the tctrahcdra 
described by these three- forms form a closed surface is then (I3.31[) . Thus the simplicial geometric reasoning 
goes through also for this new constraint. In terms of the equations of motion of the theory, on the other 
hand, the only argument for the need of this constraint is the following. If Vr" E^P = and we assume that 

Y, AB = ±e A A e B , then it follows that v!^ = as well. Integrating this equation (with the connection 
again set to zero) over a 4-ball (which can be thought of as our 4-simplex) whose boundary is a 3-sphere, 
triangulated by tetrahedra, then leads to (I3.31[) . We then however have to assume simplicity of Y> AB . A 
more direct derivation of (|3.31|) from the equations of motion would be desirable. 
The role of this constraint, anyway, is the following. Consider the closure constraint 

E AB (A AB ) + E AB (A AC ) + E AB (A AD ) + £ AB (A AE ) = 0. (3.32) 

Contracting with ua(^b) gives, using the linear simplicity constraint n J 4(A B )£^- B (A AB ) = 0, 

n A (A B )Z AB {A AC ) + n j4 (A B )E AB (A AD ) + n A (A B )Z AB (A AE ) = 0. (3.33) 

Alternatively, one may start with the 4d closure constraint and contract with S' 4b (A A b) to get, again using 
the linear simplicity constraints, 

n A (A c )^ AB {A AB ) + n A (A D )£ AS (A AB ) + n A (A E )^ AB (A AB ) - 0. (3.34) 

In total one obtains 20 + 10 = 30 equations of this kind that can be used to express some of the combinations 
n A (A)T, AB (A) in terms of others. Substituting the resulting expressions into the five discrete volume 
constraints (|3.30[) one finds that the equations (I3.30|) indeed follow from the relations (|3.33[) and (I3.34p . We 
have seen in the continuum that the summed constraints (|3.23p are sufficient to identify Y> AB = ±E A A E B , 
and hence we find that in the discrete case the situation is analogous to the case of quadratic constraints in 
that a sufficient set of volume constraints can be viewed as secondary. 

To see more clearly what happens in both our construction and in the case of quadratic constraints analyzed 
in note that in our linear case one could use the 3d and 4d closure constraints to express the variables 
n A (A E ) and E AS (A AE ), E j4B (A BE ), E j4s (A CE ), E AS (A DE ) in terms of the others. Taking the linear 
simplicity constraints into account, one is then left with twelve independent combinations n A (A)T, AB (A), 
just as in the continuum. In the continuum, we saw that one can impose the three additional constraints 
(I3.23P on the twelve contractions riAbef S^ b s to complete the identification Yi AB — ±E A A E B . In the discrete 
case, one has the following three additional conditions coming from linear cross-simplicity constraints: 

= n A (A E )E AS (A AE ) = n A (A B )S AS (A AC )+n j4 (A B )I]' 4B (A AD )+n J 4(A c )S AB (A AB ) 

+n A (A c )Y. AB (A AB ) + n A (A n )E AB (A AB ) + n j4 (A D )S AS (A AC ) (3.35) 

and similar ones coming from nA(A E )£ AB (A BE ) = and n j 4(A E )E j4B (Ac E ) = 0. These are precisely the 
analogue of the continuum constraints (|3.23|1 . 

Similarly, in the case of quadratic simplicity constraints, one can use 3d closure to eliminate 
E AB (A AE ),£ AB (A BE ),£ AB (A CE ),E AB (A DE ). Then one observes that additional quadratic cross- 
simplicity constraints give expressions such as 

= e ABC DX AB (A AE )Z CD (A BE ) = e ABCD Z AB (A AC )i: CD (A BD ) + e ABCD Z AB (A AD )Z CD (A BC \3.36) 

which are equivalent to the desired (two) volume constraints. 

All of this is an exercise in solving a system of linear equations for which there might be a more simple and 
elegant description, but the upshot is the following. The sufficient set of linear volume constraints (|3.30j) does 
indeed follow from the linear simplicity constraints and the closure constraints, once one also imposes a four- 
dimensional closure constraint on the normals to tetrahedra that seems very natural in light of their geometric 
interpretation. Just as in the formulation in terms of quadratic simplicity constraints 23[, the volume 
constraints can be viewed as secondary constraints that imply conservation of the simplicity constraints in 
time, or put differently, the volume constraints follow if the simplicity constraints hold everywhere. Once 
more this strenghtens the relationship between the discrete linear volume constraints we have identified and 
the edge simplicity constraints of [131 ] . 
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IV. LAGRANGIAN AND HAMILTONIAN FORMULATION 



Let us briefly outline the Lagrangian formulation of 4d gravity resulting from our linear constraints added 
to BF theory. One adds the linear simplicity and volume constraints to the action of BF theory using 
Lagrange multipliers: 

S = J d±x Qe^r RabM + l^eMCDVtfBSFM + H^W/E^) , (4.1) 

where the Lagrange multiplier field 'E. B bde ^ satisfies S a B hde ^ = S^ 6 ' ' de ^' , and 

e ae f'E. B bbe ^ = (no sum over b). (4-2) 
Indeed, varying with respect to then gives back the constraints 

[0, {a,b}c{d,e,f}, fA9 . 

n Adef i, ab -{ Caeff B t & = d ( forsom e/B). W 

Note that the second line corresponds to the set of constraints (|3.19|) and not to the summed version (I3.23|) . 
and that it is clearly sufficient for the geometric interpretation of Yi AB . The field equation from varying with 
respect to the connection lu is the usual 

v[: ) Sf|=0, (4.4) 

where V is the covariant derivative for the connection oj ab . The remaining equations involve the Lagrange 
multipliers, as would be expected: 

1 abcdp r ,1 ^ abed, nCDt l , ^abdef _ n ~abdef^AB _ n ( A r\ 

-e iiABcdM + -e £ABCDti c d M + =[ B riA]def — U, ^ B L ab — U. (4.5; 

We have seen that the constraints imply that Y, AB = ±E A A E B , and when substituting this back into the 
action one will recover general relativity, modulo the possible sign ambiguity we have already discussed. 

We leave a complete Hamiltonian analysis of this theory to future work. However, we note a feature of the 
theory that follows directly from the use of linear constraints, and from the introduction of the additional 
variables ua- 

As in unconstrained BF theory the initial dynamical variables will be the spatial part of the connection 
uj ab and its conjugate momentum P\ B = ^e^ k T,ABij- We also saw that the equation of motion Vr"'S^? = 
is unaffected by the constraints. Hence there will be Gauss constraints of the form 

gCD = g.pCDi + u C Et pEm + u D Et p CEi (4 g) 

on the canonical momenta. Their role is to generate G gauge transformations. 

Looking at the action (|4.ip . one would already require that (|4.6p should be modified to generate gauge 
transformations on the normals riAdef] (|4.1I) is only invariant under gauge transformations if the three- forms 
riAdef are transformed. The need for such a modification is also seen if one computes Poisson brackets 
between the linear simplicity and Gauss constraints. Define "smeared" constraints 

C[S] := J E%' def n A de f e tjk P ABk , G[A] := / A cd Gcd- (4.7) 
5C[S] SQ[A] 



One then finds that 

{m,g[A}} 



S pGHm 6cjGHm 



K ADpB D m _ A BDp Ai 



-C[A ■ S] - J ~f def A AD n Ad e f e t]m P B D m , (4.8) 
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where (A-S)^' de ^ = Ajj B E l ^' de ^ . The first term alone would imply that G[A] generates gauge transformations, 
but the second term is an unwanted extra piece. For Q to be a generator of gauge transformations, it must 
be first class (i.e. commute with other constraints up to linear combinations of constraints). We can remedy 
this by adding the variables riAabc to the phase space, together with their conjugate momenta n Aabc . This 
extension of the phase space is analogous to the situation considered in |18l sect. 3] for a generalised quadratic 
constraint formulation of Plebanski theory. Now we can define a new Gauss constraint 

g,CD = gCD _ ^Dlate (4 J) 

Then, computing the Poisson brackets of the new Gauss constraint with C[H], one finds 

{C[H],S'[A]} = {C[S],£?[A]} - J ^ def A CA n Cdef e ijk P ABk = -C[A ■ S], (4.10) 
as desired. We have however increased the number of phase space variables at each point by 32. 



A similar reformulation of the Gauss constraint, leading to a relaxation of the gauge invariance properties 
of spin network states, has been already suggested by the Hamiltonian analysis of the Plebanski theory 16[ , 



and it has been advocated in the loop quantum gravity context in [39, 
field theory context 0, S, |H| . 



40j as well as spin foam and group 



V. SUMMARY AND OUTLOOK 



We have investigated a formulation of classical BF-Plebanski theory where the constraint T, AB = ±e A Ae B , 
needed to reproduce general relativity in four dimensions, starting from topological BF theory, is imposed 
through constraints linear in the bivector field S AB . The discrete counterpart of a part of these linear 
constraints (th e 'simplicity constraints'), in fact, has proven very useful in the spin foam approach to quantum 
gravity 



21, 22, 2 



The corresponding continuum constraints have been easily identified, and can indeed be used to replace 
the quadratic "diagonal" and "off-diagonal" parts of the simplicity constraints appearing in the Plebanski 
formulation. As in the discrete case, one needed to introduce a new set of variables ha which are assumed 
to form a basis of three-forms at each point of spacetime, and are slightly stronger than the quadratic 
constraints: they eliminate two of the four sectors of solutions that are present for quadratic constraints. 

In the second part of the analysis we found that the quadratic volume constraints of the Plebanski formu- 
lation, needed to complete the identification T> AB = ±e A A e B , can also be replaced by linear constraints, 
which again are stronger than their quadratic analogues. They do not require an additional non-degeneracy 
assumption on T, AB . However, a non-degeneracy assumption on the three- forms tia is still necessary, and 
only when this is imposed one can hope to eliminate all "non-geometric" degenerate configurations for Yj AB , 
which are feared to dominate the quantum theory in the case of quadratic volume constraints. Also, while 
for quadratic volume constraints the variables tia merely specify normals to submanifolds {x a — constant} 
and hence can be independently rescaled arbitrarily at each point, for linear volume constraints they directly 
specify, up to an overall rescaling, the frame field encoding the metric geometry, i.e. an orthonormal basis 
in the cotangent space at each spacetime point. 

We have then analyzed the discrete (simplicial) translation of the linear constraints we identified. In 
the context of spin foams, the quadratic volume constraints follow from imposing the (quadratic) diagonal 
and off-diagonal simplicity constraints everywhere together with closure constraints on the discrete variables 
Yi AB (/\). We have shown a similar property for the linear volume constraints. If (linear) diagonal and 
off-diagonal simplicity constraints and closure constraints for both bivector variables T, AB (A) and normals 
n^(A) are imposed everywhere, a sufficient set of linear volume constraints follows. This means that "non- 
geometric" bivector configurations cannot appear if the additional closure constraint on the normals holds, 
and the same normals are assumed to be non-degenerate. 
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We have not performed a complete Hamiltonian analysis of the resulting linear constrained BF action 
for gravity, but only noted that the use of linear simplicity and volume constraints immediately requires 
a modification of the usual Gauss constraint to generate a transformation of the normal 3-form variables 
rbA alongside that of the E's; a similar relaxation of the Gauss constraint, which translates at the spin 
foam and discrete gravity level into a closure constraint for simplices, and in the canonical quantum gravity 
context into a generalisation of spin network states, has been suggested on more than one occasion in the 



literature 26, 39 41], even if its proper implementation at the quantum level has not been yet developed. On 



the classical level, therefore, a full Hamiltonian analysis of the constraints would be highly desirable. This 
would involve adding momenta for the components which are Lagrange multipliers in unconstrained 

BF theory, as well as those for the normals tia we have introduced, so that all variables can transform 
nontrivially under G gauge transformations generated by a modified Gauss constraint, as shown. 

Still at the classical level, but with obvious implications for the quantisation, one aspect of our construction 
that deserves further work is the relation between the discretised linear volume constraints we have found 
and the edge simplicity constraints used in 13] , in turn related to the conditions on connection variables of 



37] . As noted, the two sets of constraints appear to be very similar, and their role in the classical theory is 



the same, in particular, they remove (partly) the non-geometric configurations from the configuration space 
(or phase space) of the theory and appear as "secondary" in the sense specified above. So it natural to 
conjecture that one is simply a reformulation of the other. The implications for the quantum theory are not 
only due to the dominant role that non-geometric configurations may play in the quantum theory, if not 
removed, but also in the fact that one discrete formulation of these constraints can actually be simpler to 
implement in a spin foam context than the other. 

The possible use of our findings in the spin foam and group field theory context, and more generally in any 
quantisation based on the formulation of gravity as a constrained BF theory, are in fact most interesting. In 
particular, it seems to be important to explore how a closure constraint on normals could be implemented 
into existing spin foam models, given that we found it to be necessary for the full imposition of the geometric 
constraints on the variables of topological BF. A convenient setting to do so could be the GFT formulation 
of (2^ |. since there the simplicial geometry and the contact with classical actions is brought to the forefront. 
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